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OO ■ Abstract 

o: 

. We study the collapse of a spherically symmetric dust distribution in d-dimensional AdS spacetime. 

' We investigate the role of dimensionality, and the presence of a negative cosmological constant, in 

O . determining the formation of trapped surfaces and the end state of gravitational collapse. We obtain 

D ' the self-similar solution for the case of zero cosmological constant, and show that one cannot construct 

a self-similar solution when a cosmological constant is included. 

\o ■ 

1 Introduction 

' There are many models of spherical gravitational collapse in classical general relativity which exhibit the 
f^. formation of black holes as well as naked singularities, starting from regular initial data [I], [2]- The study 
I of quantum effects in the vicinity of the gravitational singularity then becomes significant. Such studies can 
CN ' be divided into two classes : (i) quantum field theory in curved space, and (ii) quantum general relativistic 
treatment of gravitational collapse. 

The earliest investigations of quantum field theory in the dynamical background of a collapsing spherical 
star were probably those due to Ford and Parker [3] and Hiscock et al. [1]. These works introduced 
important techniques, such as the calculation of the quantum flux in the geometric optics approximation, 
, and the regularization of the 2-d quantum stress tensor, which were used extensively in later studies. A 
I systematic study of semiclassical effects in gravitational collapse was initiated by Vaz and Witten in [5], 
O ■ [6] and pursued in a series of papers [7], [8], [9], [10], [11], [12], [13]. Typically, these studies showed an 
J> important and interesting difference in the nature of quantum particle creation between the two cases - one 
in which collapse ends in a black hole, and another in which it ends in a naked singularity. The formation 
' of a black hole is accompanied by the emission of Hawking radiation, as expected. However, when the 
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collapse ends in a (shell-focusing) naked singularity, there is no evidence of some universal behaviour in 
the nature of quantum emission. It is typically found though, that the emitted quantum flux diverges in 
the approach to the Cauchy horizon. This divergence disappears when the calculation of the quantum 
flux is terminated about a Planck time before the formation of the Cauchy horizon, when the semiclassical 
approximation breaks down. Instead of the divergence, one finds that only about a Planck unit of energy 
is emitted during the semiclassical phase, and a full quantum gravitational treatment of the physics of the 
singularity and the Cauchy horizon becomes unavoidable. These developments have been reviewed in |14j . 

A full quantum gravitational treatment of collapse can be performed via a midisuperspace quantization 
within the framework of quantum general relativity. The aims of such a programme are manifold - to 
construct a quantum gravitational description of the black hole; to check if the gravitational singularity can 
be avoided in quantum gravity; to obtain a statistical derivation of the black-hole entropy from quantum 
gravitational microstates; and to determine the role of quantum gravity in ascertaining the nature of 
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quantum emission from a naked singularity. The midisuperspace quantization programme has been carried 
out by us in a series of papers [15] , [16] , [T7] , [18] , [19] , [20] , [21] , and work along these lines is still in progress. 
It is fair to say that while some progress has been made on aspects related to quantum black holes and black 
hole entropy, issues related to singularity avoidance and the nature of quantized naked singularities have 
thus far proved difficult to address, largely bcause of problems relating to finding a suitable regularization 
scheme for the quantized Hamiltonian constraint in canonical general relativity. Also, we still do not have 
a definitive answer as to the nature of quantum gravitational corrections to the semiclassical spectrum of 
Hawking radiation. By this we mean the following: starting from a candidate theory of quantum gravity 
such as quantum general relativity, one can derive Hawking radiation in the semiclassical approximation 
using a suitable midisuperspace model. Going beyond the semiclassical approximation, it is expected that 
quantum gravity will induce (possibly non-thermal) corrections to Hawking spectrum, but this still remains 
to be worked out in its full generality. It is hoped though that some progress will be possible on these 
unsolved problems if one makes contact with the methods of loop quantum gravity. 

All the classical and quantum studies mentioned so far have pertained to gravitational collapse in 3+1 
dimensions. Motivated by the desire to overcome some of the obtacles faced in 3+1 physics, we turned 
attention to investigation of 2+1 gravitational collapse. Homogeneous dust collapse in 2+1 dimensions 
was first studied in [22] and for the case of collapsing shells in [23J . This lower dimensional model, though 
simpler in some aspects, throws up new fascinating issues of its own, which have been studied in the 
context of inhomogeneous dust collapse in [25], [2^, [27]. Classical 2+1 collapse admits a naked singularity 
for some initial data, but there is no coresponding quantum particle creation. A black hole solution (the 
well-known BTZ black hole [28]) is possible in the presence of a negative cosmological constant, but the 
thermodynamics and statistics of the quantized BTZ black hole is completely different from that of the 
4-d Schwarzschild black hole. 

These differences prompt us to the following question: in determining the nature of thermodynamics 
and statistics of the quantized black hole, and the nature of quantum emission from naked singularities, 
what is the role of the cosmological constant, and of the number of spatial dimensions? The present paper 
is the first in a series of three papers which addresses this question, by studying classical and quantum 
aspects of spherical dust collapse in an AdS spacetime with an arbitrary number of dimensions. In the 
current paper, we solve the Einstein equations for a collapsing dust ball in an asymptoticaly AdS spacetime, 
and examine the nature of the gravitational singularity. Quantization of this model will be taken up in 
two subsequent papers. 

The plan of the paper is as follows. In Section II we give results for spherical gravitational collapse 
of dust in an asymptotically flat d-dimensional spacetime. While this problem has been studied earlier 
by various authors [29] , [30] , [31] , [32] , we present here a simpler derivation of the occurrence of a locally 
naked singularity, and also obtain new results on the self-similar solution. More importantly, the results of 
this Section serve as a prelude to the corresponding analysis presented in Section HI, for collapse in an AdS 
spacetime with arbitrary number of dimensions. While gravitational collapse of dust in four dimensional 
spacetime with a positive cosmlogical constant has been studied in [33], (see also [33]), and for a negative 
cosmological constant in |24j , to the best of our knowledge dust collapse in a d-dim AdS spacetime has not 
been studied before. 

One could question the introduction of a negative cosmological constant, as is done in this paper, when 
the observed Universe has a cosmological constant which is perhaps positive, or at best zero, but certainly 
not negative. Firstly, collapse physics in a deSitter spacetime is complicated by the presence of a deSitter 
event horizon, in addition to the black-hole event horizon. It thus seems natural to first address the AdS 
case before moving on to the more realistic, and more difficult, deSitter case. There are also reasons 
to believe that it would not make sense to directly construct a quantum black hole model in a higher 
dimensional space with a positive cosmological constant, because quantum gravity in such a spacetime 
may not exist nonperturbatively [35], [36] . Pure quantum gravity with a positive cosmological constant 
may hence not exist as an exact theory, but only as a part of a larger system [36] . It is also a question of 
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great interest as to whether studies of statistical properties of AdS black holes in canonical quantum general 
relativity can benefit from what is known about the AdS/CFT correspondence, as suggested recently in 
[37j for the 4-d case. 



2 Higher Dimensional Spherically Symmetric Dust Collapse in the Ab- 
sence of a Cosmological Constant 

2.1 Solution 

The metric for a spherically symmetric space-time can be written in the form 



where 



-del + sin^ 9i{dei + sin^ ^2(d^i + + sin^ On^idel)). 



(1) 



(2) 



Here the number of space-time dimensions is (n+2) where n > 1 is the number of angular coordinates 
and the 2 designates one time dimension and one radial dimension. For the case where the cosmological 
constant A = 0, Einstein equations are 

Gfiu — kT^y (3) 

where is a constant related to Newton's constant of gravitation G (see section (2.4)) and T^y is the stress- 
energy tensor. For the case of non-rotating dust one can choose a synchronous and co-moving coordinate 
system in which the only non-zero component of the stress-energy tensor is Too = ^it^T)-, where €{t,r) is 
the energy density of the dust. Further, in co-moving coordinates the goo component of the metric can 
be chosen to be minus one. With this choice for the metric in ([T]) we get the following independent set of 
Einstein equations 
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Components G33, G44 etc. are given by expressions similar to that for G22 except for overall sine squared 
factor (s). The Ricci scalar is given by 

n = ^[-2n{n-l)(^R'^ + e^{l + R^)^ +2nR(^R'X' -2R" + e^{RX + 2R)^ 

+e^R\X^ + 2X)]. (8) 



Solving the equation for Goi we obtain 

l + /(r 



= (9) 
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R' = m + l^v (10) 



In the above expression /(r) is an arbitrary function called the energy function. Integration of the equation 
for Gil after using equation ([9]) gives 

F(r) 

Here F{r) is another arbitrary function and is called the mass function. In what follows we will only 
consider the so called marginally bound case for which /(r) = 0. In this case (llOp can be integrated easily 
and after choosing the negative sign for the square root corresponding to in-falling matter we get 

n + l 

t-tcir) = ^-^^ (11) 

where tc{r) is yet another arbitrary function which can be fixed by using the freedom in the choice of the 
r-coordinate. We relabel r such that at t = 0, i? = r. With this choice we have 

n + l 

Ur) = — (12) 

From the above equations we see that at t = tc(r), R{t,r) = and this implies singularity formation for 
the shell labeled r as indicated by the blowing up of the Ricci scalar in ([5]) . Finally, substituting for A from 
([9]) in the equation for Gqq we find that 

n F' 

ke{t,r) = (13) 
From this one can obtain an expression for the mass function 

F(r) = — / e(0,r)r"dr. (14) 

2.2 A simple derivation of the naked singularity 

We now look at the nature of the R = singularity formed at the center r = of the dust cloud. For 
this we follow the method used in [38] and start by assuming that the initial density profile e(0, r) has the 
following series expansion near the center r = of the dust cloud 

e(r) = eo + eir + — r + .... (15) 

Using this in ()14p we find that in this case the mass-function can be written as 

F(r) = F„+ir-+i + F^+ar'^"^' + i^n+3r"+=^ + (16) 

where it is to be noted that n is not a free index but, as before, refers to the number of angular dimensions 
and 

F - ^'-^ h7\ 

^^+^- n{n + i){i -1)1 ^ 

and i = 1, 2, 3.... From (I12p we know that the singularity curve is given by 

n + l 

2 r— 

t,{r) = ^ (18) 

The central singularity at r = forms at the time 



° n + l^TWr n + lV 2/eeo ' ^ ^ 
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Here, as a special case, we note that when e(r) = eo, a constant (Oppenheimer-Snyder collapse), F{r) = 
Fn+ir^'^^ and the singularity curve is given by ts = 2/(n + l)y/Fn+i which is independent of r implying 
that all shells become singular at the same time as the central shell. Near r = one can use the expansion 
for F(r) as in (|16p and approximate the singularity curve as 
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In the above equation F^+i is the first non- vanishing term beyond -Fn+i in the expansion for F{r). 

One would like to know whether the singularity at i = toi = is naked or not, and for this we focus 
attention on radial null geodesies. We want to check if there are any outgoing radial null geodesies which 
terminate on the central singularity in the past. Assuming that there exist such geodesies we assume their 
form near r = to be 

t = to + ar"" (21) 

Comparing this with (|20p we conclude that for the null geodesic to lie in the spacetime one must have 
a > i — 1 and if a = i — 1 then 

Fn+i 



a < 



(22) 



(This is because Fn+i is negative, which will be the case if we demand that e(0, r) be a decreasing function 
of r). 

Since one is interested in the region close to r = 0, we expand (jlip to leading order in r to obtain 
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From the metric one finds that for null geodesies dt/dr\i\jQ = R' . Diff'erentiating (|23p w.r.t. r we get 
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Along the assumed geodesic, t is given by ([211) . Substituting this in (p^ and equating it with the derivative 
of (j2ip . i.e. dt/dr = aar°^~^ gives 



aar 



a-l 



[1 
[1 



(n + l) 
2 

(n+l) 



n+l 



IF 



2 Fn+1 

/FU7(to + ar") - ■ 



("-!) 



(to + ar°)] 



n + 2i - 1) Fn+i i-i 



n+l 



(25) 



This is the main equation. If it admits a self-consistent solution then the singularity will be naked otherwise 
not. To simplify this we note that \/Fn+ito = 2/(n + 1), as follows from (119p . 
We first consider the case a > i — 1. To leading order this gives 
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This equation implies 

n + 1 ' V 2F„+i; ^ ^ 

Since .Fn+i is the first non-vanishing term beyond i^n+i; we have the condition i > 1. Also, for consistency 
we require a = {n + 2i — l)/{n + 1) > i — 1, which together with the previous condition on i implies 

1 <i < 2n/{n - 1). 

This implies that in 4-dimensions, where n = 2, we have 1 < i < 4, which means that i = 2,3 are 
the allowed values. That is, models for which either ei < (corresponding to i = 2) or ei = 0, €2 < 
(corresponding to i = 3) will have a naked singularity. 

Similarly in 5-dimensions, where n = 3, we find that 1 < z < 3 implying i = 2, i.e. only for ei < 
we get naked singularity. In 6-dimensions, n = 4 and we have 1 < i < 8/3, implying i = 2 as the 
only allowed value, i.e. the singularity is naked only if ei < 0. One notes that for all higher dimensions 

2 < 2n/(n — 1) < 3 and therefore only i = 2, i.e. ei < gives naked singularity. 

As another special case we note that for n = 1, that is in (2+1) dimensions, a = (n + 2i — l)/(n + 1) = i 
and therefore the condition a > i — 1 is always satisfied, implying that in this case we always have a naked 
singularity, which is in agreement with what has been observed in earlier work on (2+1) dimensional dust 
collapse |25j. 

We next consider the case where a = i — 1. Here (f25]l gives 

(n-l) 

(i-l)ar = — VFn+ia 



2 V 2F, 



■t+i 



(n+1) r- (n + 2z - 1) F„+A 

which implies i = 2n/{n — 1). Now the conditions on i are that it be an integer greater than 1. These 
two conditions are met only for n = 2 and n = 3, that is, in (3+1) dimensions and in (4+1) dimensions 
respectively. For n = 2, i = 4 (which corresponds to £3 < 0) and for n = 3, i = 3 (corresponding to £2 < 0). 
Since the 4-dimensional case, corresponding to n = 2, is already reported in the literature [39] we focus 
attention on the 5-dimensional case corresponding to n = 3. 
Substituting n = 3, i = 3 in (p8]) we obtain 

8v^a3 + + 4F4) + ^« + f = 0- (29) 

The above cubic for a has to be solved subject to the constraint < a < —F^/AFJ as mentioned earlier. 
By defining o = y/F^b and F^ = F^^ the above equation is simplified to 

26^(46 + + (2fo + 0^ = (30) 

and the constraint on a results in a constraint on b given by < 6 < — 1^/4. By defining —h/^ = Y and 
—l/^ = r], the above cubic is further simplified to 

2Y'^{AY -l)+r]{2Y -if = 0. (31) 

For a naked singularity to form this equation for Y should have a positive root subject to the constraint 
< y < 1/4. 

Now for a general cubic 

qqx^ + 3aix^ + 3022; + 03 = 0, (32) 
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if we define H = aoa2 — a\ and G = CqOs — 3aoaia2 + 2af , we have the following conditions on the roots 
of the cubic [40J : 

1. + 4:H^ < 0, the roots of the cubic are all real. 

2. + AH^ > 0, the cubic has two imaginary roots. 

3. G^ + AH^ = 0, two roots of the cubic are equal. 

4. G = and H = 0, all three roots of the cubic are equal. 

Using these we can find the conditions on r/ for which the cubic in (jSip has at least one real root in 
the desired range. Here it should be noted that r/ as defined above has to be positive. It is found that for 
< r/ < (—11 + 5\/5)/4 all the three roots are real and at least one of these satisfies < y < 1/4. For 
r] > (—11 + 5\/5)/4 the real root is negative. The range of r/ found above implies that for ^ < 4/ (11 — 5\/5) 
one gets a naked singularity. 

We also note that for the Oppenheimer-Snyder collapse mentioned earlier, no naked singularity is 
formed since all shells become singular at the same time. 



2.3 Formation of Trapped Surfaces 

We now consider the formation of trapped surfaces. For this consider a congruence of outgoing radial null 
geodesies with tangent vector = dx^/dk where /c is a parameter along the geodesic and i = (0, 1) |41j . 
The expansion for these geodesies is given by 



From this one finds that 
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Trapping occurs when 9 = and the above equation with i?' > 0, i? > and K"^ > implies that this 
condition is met for 

R^ = ^- 

In 4-dimensions where n = 2 we get the well known result 



^ = 1 (36) 

For the general case one finds that the time at which trapping occurs ttr is given by 

2 / r~ 1 \ 

ttrir) = -— ^= - F{r)—r (37) 



n 



+ 1 



which means that the central shell is trapped at ttr{r) = 2/(n + l)-^F„+i, that is, at the same time as 
the formation of the central singularity. For the outer shells trapping occurs before those shells become 
singular. 



2.4 Exterior Solution and Matching with the Interior 

We take the metric in the exterior to be independent of time and given by 

ds"^ = -f{x)dT^ + g{x)dx^ + x^dQ^ (38) 
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where {T,x,9i,62---) are the coordinates in the spacetime exterior to the dust cloud. The components of 
the Einstein Tensor corresponding to the above metric are 
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Gil 
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and G33, G44 etc. are the related to G22 as in the interior. Solving the vacuum Einstein equations G^^i, = 
one finds 



g{x) = 1 
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Here C is a constant of integration. Thus the exterior metric is the Schwarzschild metric 
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For this to be a valid solution in the exterior we need to match the metric coefficients as well as their first 
derivatives (extrinsic curvature) in the exterior with the corresponding quantities in the interior at the 
boundary of the dust cloud r = Vg say, [42] , [43] . This will also determine the only unknown quantity G in 
the Schwarzchild solution. At the surface the exterior coordinates will be some functions x = x{t, r^) = Xs{t) 
and T = r(t, r^) = Ts{t) of the interior coordinates. These relations imply dT = Tgdt and dx = Xsdt. 
Therefore at the surface (where dr = 0) 
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dt^ + xld^f = -dt^ + Rl{t)dn 



(45) 



Matching the metric coefficients for gives Xs{t) = Rs{t) and matching the metric coefficients for dt^ 
then implies 



G 
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1. 



(46) 



To match the extrinsic curvature (second fundamental form) we need the normal to the surface. In the 
interior coordinates the components of the normal are found to be = (0, -R', 0....0). Similarly in the 
exterior coordinates the normal is given by = (— i?^, T^, 0...0), where the relation dx — Rgdt = was used. 
The extrinsic curvature is given by Kab = nfj_-^eae'^, where Ca = dx^/dy'^ with x^ being the coordinates of 
the (n+2)-dimensional manifold and y"" being the coordinates on the boundary of the manifold. Since there 
is only one undetermined constant G, we match only the Kq^q^ component of the extrinsic curvature. It 
can be easily checked that the other components do not give anything new. We find that at the surface the 
extrinsic curvature in the interior coordinates is given by -f^g^e^ = Rs- Similarly in the exterior coordinates 
we have Kg_^g_^ = i?^ (l — G/R'^^^) Tg. Equating these two expressions for Kq-^q-^ gives 
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Using (P6j) and Rg = F{rs) / R^~^ (see (fTOj) ) the above equation gives C = F{rs), where from (fH|) it is 
clear that -F(r<j) is proportional to the total mass of the dust cloud. Thus we find that for the metric 
coefficients and their first derivatives to be continuous across the boundary the metric in the exterior is 
given by 

ds^ = -{l-^) dT- + (l - " + -'dn- (48) 

Now F{rs) = (2A;/n)/o' e(0,r)r"dr and we know that mass of the dust cloud is given by M = 
/q" e{Q,r)dV where dV is the volume element of a spherical shell lying between r and r + dr \n {n + 1) 
space dimensions. This volume element is given by 

n + l 

dV = -^r"dr (49) 

T^f n+l \ \ ' 



Therefore 

n + l 

27r— 



M = -—^ e{0,r)r-dr. (50) 
i [ — ) Jo 

e{0,r)r-dr = (51) 



This implies 

/o ■ ■ 2tt"^ 
Using this we find that the mass function can be written as 

" 27r— 

One can also find the constant C in the Schwarzchild solution using the weak field limit. For this 
we assume that Newton's law for gravity holds for any number of dimensions i.e. V ■ g = 47rGe(r) and 
g = — V(/)(r). Here g is the gravitational field strength and (j) is the gravitational potential (note: Newton's 
gravitational constant G being dimensionful will be different in different dimensions, however, this does 
not affect the form of the equations). Using this we find that in (n + l) spatial dimensions the gravitational 
potential is given by 

2{n - l)7r — r"-i 

where n > 1 (potential has a logarithmic dependence on r in 2 + 1 dimensions). In the weak field limit 
the Schwarzchild solution is goo ^ — (l — and gu — > (l + :^:^Lt) ■ Also, using geodesic equation we 

find that generically, in the weak-static field limit goo = ~(1 ~ 20) and gu = (1 + 2<j)). Comparing the two 
expressions for (700 (or for gu) one finds that C = 20r"~^ and using the expression for cj){r) as found above 
one gets 

(n — 1)TT 2 

This expression for C will be the same as that found above from matching if the constant in Einstein's 
equations is chosen to be A; = 4mTG/{n — 1). For n = 2 this reduces to the value SvrG as used in 
4-dimensional theory and which when used in (j48p results in the familiar Schwarzchild solution 

= _ (1 _ ?|11) + (1 _ ?|1£) <i,^ + (55) 
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2.5 The Self-Similar Solution 

To see the effect of dimensions on the nature of quantum particle flux (which will be described in a work 
subsequent to this), we would like to have a globally naked singularity. It is known that a locally naked self- 
similar solution is also globally naked }44j . where self-similar spacetimes are defined by the existence of a 
homothetic Killing vector field. Therefore here we look at the dependence on dimensions of the self-similar 
dust model. 

In a self-similar collapse any dimensionless quantity made from the metric functions, has to be a function 
only of t/r. This can be seen by starting from the definition of a homothetic Killing vector field ^, 

Ca;/3 + = 25fa/3 (56) 

(We emphasize here that we are dealing only with self-similarity of the first kind, which is defined by the 
above equation). This condition implies that 

C-iG^, = 0. (57) 



where G^^ is the Einstein tensor for the spacetime [45] , [46] . It follows then that the the energy momentum 
tensor should also satisfy the equation 

C(T^u = 0. (58) 

For a perfect fluid the energy-momentum tensor is T^y = {p + lj)u^Uy +pgf_iu- The condition (j58p implies 
the following, 

C^u^' = -u^ (59) 

= -2^ (60) 

L^p = -2p (61) 

Based on the above equations, the existence of a homothetic vector implies that the metric components in 
a comoving coordinate system will be of the form such that the dimensionless quantities become functions 
of r/i, [U]. Following [¥7] we assume a general spherically symmetric ansatz of the form 

ds^ = -e^'^dt^ + e2*dr2 + R^dQ^ (62) 

and also we have 

C^u^' = -u^ (63) 

If we assume the vector has only r and t components given by = abr + ^"^^ expand (I56p we get 
PDEs for the vector components. The condition of comoving metric implies the condition that a^t = 

and f3^r = 0. After redefining the independent variables to = and i^t = t/f3 we now redefine the 
dependent variables by 

§ = ^ + log(a) - log(f) (64) 

= + log (/?)- log (f). (65) 
Under the above change of variables the equations (f56]) become. 



rRf + tRi = R (66) 

f^f + i^t = (67) 

r^f + t^t = (68) 

So this shows that we can go to a coordinate system in which the metric functions are functions of f/t . 
So ^ and $ are functions of z and R is rTZ{z). 
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It can now be shown that for spherically symmetric, self-similar dust collapse, the mass function is 
given by -F = Ar"~^, where A is a constant. For this we start with the Gqi component of Einstein equation 



G, 



01 



n {R'\ - 2R') 
2 R 



0. 



(69) 



Defining the self-similarity parameter z = r/t and writing R = rR, where R and A being dimensionless are 
functions only of z, the above equation can be written as 



^dX ^dR ^dXdR ^_^dPR_^^ 
dz dz dz dz dz"^ 



(70) 



where we have used R' = R + zdR/dz and R = —z^dR/dz. The above equation is solved easily to obtain 



c{R + z 



dz 



(71) 



where c is a constant of integration and equals one for the marginally bound case (and will therefore be 
ignored in what follows). Similarly the Gu component of the Einstein equation gives 



n{n-l) ~ dR\ nin - I) ~ dR\ 

— R + z— ^ R + z— 

2 \ dz 2 \ dz 



1 + z^ 




-nzR \ R + z 



Solving this we obtain 



ir-^R!^ 



where c is a constant of integration. 

Now the mass function is given hy F = R^R^^^ and using R = —z^dR/dz this can be written as 



Using ()73p in the above equation we obtain F = cr^ ^ which is the desired result. 
With this result (jlip becomes 

1+1 n + 1 



R— 



VAr"-i(6'r - t) 



dR\ / dR ^3 d^R 
dz I \ dz dz"^ 



(72) 
(73) 

(74) 
(75) 



Here 6 = tc/r is a constant. This is because in self-similar collapse any dimensionless quantity has to be 
a function only of t/r whereas tc/r, being a function only of r (see (|12|) ) has to be a constant. We are 
interested in finding the behavior of density e = nF' /2kR' R"" in the neighborhood of the centre r = 0. 
Using (|75|) and F' = A(n — l)r"'~^ in the expression for density we find 



4n 



n — 1 



2n 
n — 1 



etr + f 



K{n + 1) 

For r ~ we neglect the second order term in the above equation and obtain 

n -1 



4n 



K(n+ l)t2 



1 



2n r 



n-lt 



(76) 



(77) 



Also for r ^ 0, R^~^^ = Xt^r^ which implies 



{n + l)t\^ 



n + l 

R~. 



(78) 
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Substituting this in ([77|) we get 



4n 



1 + 



2ne 



n 



1 



2 

an"' + 



n + l 



n 



1 



' 4 2(n + l) 



(79) 



where a„ = 2/(n + 1)\/A and z = R/t. This shows how the density profile should depend on the number of 
dimensions to obtain a self-similar solution. The above form for density profile implies that in 4-dimensions 
(n = 2) and in 5-dimensions (n = 3) the self-similar solution corresponds to an analytic density profile 
whereas in higher dimensions the density profile is no longer analytic. 

Here we also note that for n = 1, that is in 2 -|- 1 dimensions, F = X and is thus independent of r 
and therefore one requires that energy density e should be zero. Thus self-similarity in 2 -|- 1 dimensions is 
inconsistent with the presence of matter. 



3 Spherically symmetric inhomogeneous dust collapse in the presence 
of a negative cosmological constant 

3.1 Solution 

In the presence of a cosmological constant A, Einstein equations are given by 

G^^ + Agf,u = kT^^. (80) 

For the case A < we take A —A in which case the Einstein equations become Gf_iu — Ag^y = kT^y^ 
where now A > 0. The expressions for the components of Einstein tensor are still the same as in the A = 
case. In particular since qqi = 0, therefore we again have 

a,,^^^!^^o. (81) 

The solution of this equation is again given by ^ and we again consider the marginally bound case so 
that /(r) = 0. The 1-1 component of Einstein equations is 



1 



n(n — 1) / ,9 5, , ■ r,\ 5, 



2 

Integration of this equation gives 



R'-e^{l + R^)] -ne^RR 



2A o F(r) 



AR!"^ = 0. (82) 



i?2 = + 1^ (83) 

n{n+l) ^ ' 

where as before F{r) is the mass function. Integrating this equation after taking the negative sign for the 
square root (to account for in-falling matter) we get 



2sin-i I^^B-±1^ 



, N V n(n+l) F . 

t-tc{r) = ' ■ (84) 



Relabeling the r coordinate as in the previous case so that at t = 0, i? = r we get 



2sin-i / 2A r"+i 



Ur) = ■ (85) 

For t = tc{r) we again get R{t, r) = corresponding to the singularity formation for shell labeled r. From 
Goo — ^500 = ke{t,r) we obtain an expression for F{r) which is again given by p^ . 
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3.2 A simple derivation of the naked singularity 

As before we want to see if null geodesies can come out of the singularity. For this we proceed as before 
assuming that the density profile near the center is given by 



e(r) = eo + eir + e2-^ + ... 

From the form of the mass function F = ^ J e(0, r)r'^dr we have 

F(r) = F„+ir"+i+F„+2r"+2 + .... 

where Fn+i = ni^+i) ' -^'^o™ ([8^ and ([85]) we see that the singularity curve is given by 



(86) 



(87) 



sm 



tJr) 



2A r"+i 
n(n+l) F{r) 



n+l 



2A 



(88) 



i(n+l) 



This implies that the central singularity at r = forms at time 

to -- 



2n 



2A 



(n + l)A 



sm 



n{n + l)F„+i 



(89) 



We again note that, as in the A = case, when e is a constant all shells become singular at the same time 
as the central shell. 

We now rewrite the expression for the singularity curve as 



sm 



(n + 1)A 
2n 



ts{r) 



2A 



n{n + 1) F{r) 



(90) 



It is reasonable to assume that for shells near r = the time for singularity formation is close to the time 
for the central shell to become singular i.e. ts{r) « to ^ind we can therefore write ts[r) = Ats(r) + to where 
because of the assumption made /\tg{r) ~ 0. Using this we expand the left hand side of the above equation 
using the addition formula for sines and make use of lim^j^o sin(x) = x and \\m.x^o cos{x) = 1 to get 



:ra + l )A 
2n 



cos 



[n + l)K 
2n 



to 



sm 



(n + 1)A 
2n 



to 



F 



2F, 



n+l 



(91) 



Here F^+i is the first non-zero term beyond Fn+i and is negative since we assume a decreasing density 
profile. Using Ats{r) = ts{r) — to in the above equation we can finally write the expression for singularity 
curve for shells near the center as 



ts{r) = to 



2n 



(n + l)A 



tan 



(n + l)A 
2n 



to 



F 



2F 



n+i ^i-l^ 
ra+1 



(92) 



To know whether the central singularity at t = to, r = is naked or not we focus attention on radial 
null geodesies and check if there are any outgoing radial null geodesies which terminate on the central 
singularity in the past. We proceed as in the earlier case, assuming that there exist such geodesies and 
take their form near r = to be 

t = to + ar"^ (93) 
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where, comparing with (j92p . we see that a > i — 1 and if a = i — 1 then 

2n l(n + l)A_ F^+i 



a < 



(n + 1)A 



tan 



(n + 1)A 
2n 



to 



2K 



n+l 



(94) 



for the assumed geodesic to he in the spacetime. We use (I85p and ()87p (retaining only the first two non-zero 
terms in the latter in the r ~ approximation) in (j84p to get 



n + l 



n(n + 1) 
2A 



F, 



n+l 



sm 



2n 



2n 



n+l 



(95) 



Near r = 0, the time t appearing in the geodesic equation satisfies t k, and therefore the argument of 
the sine function in (|95p is close to zero and we use the approximation sinx ~ x obtaining 



R 



Xn + l) 
^ 2 

(n + l) 



/-= / nin + 1) 

VF^+ito - \l — — tan 



2A 

n+l) Fn 



:n + l)A^J^^.„, 



2n 



2^ 



n+l 



--tor' 



I (n+l) Fn+. ^.-l^j^ 



+1 



n+l 



(96) 



From the form of the metric we know that the radial null geodesies satisfy dt/dr\MQ = R' . We take the 
spatial derivative of the above equation, substitute for t from (I93p and equate the result to the derivative 
of 



aar 



a-l 



n + 2.-l)./n(n + l)^^JJ(n + l)A^^ 



(n + l) V 2A 
(n + 2i - 1) Fn 



2n 



(n + l) 



2\fF^+i 



n(n + 1) 



2A 



tan 



:n + l)A 
2n 



to 



F, 



2\fF\i 



n+l ^i-1 



+1 



(n + l) 



-lar 



2 ' ' 2^Fn+l 

Consider first the case a > i — 1. Keeping terms only to lowest order in r we get 



aar 



a-l 



From this we have 



a 



n + 2i-l 
n + l 



n + 2i - 1 
n + l ■ 



"(" + ^)tanf,/(" + ^)\ 



2A 



F 



n+i 



2n 



2VK 



+1 



2 

n+l 



■2(i-l) 



F„. 



n{n + 1) 



2 

n + l 



2 VTUT V "(^ + - 2A 



(97) 



(98) 



(99) 



where we have substituted for in the argument of tan. Since the form of a is exactly the same as in 
the A = case we find that the conditions for naked singularity formation are also the same as mentioned 
after (j26p . That is, in 4-dimensions (n = 2), 1 < i < 4 implying i = 2,3 are the allowed values so that we 
get naked singularity for ei < or for ei = 0, £2 < 0. Similarly in 5-dimensions (n = 3), 1 < i < 3 implying 



14 



that only i = 2 is allowed so that we get naked singularity only for ei < 0. In all higher dimensions we get 
naked singularity only if ei < 0. 

Again n = 2, i = 4 and n = 3, i = 3 are critical cases satisfying a = i — 1. To analyze these we proceed 
as in the A = case. For a = i — 1 (1971) becomes 



(i — l)ar 



i-2 



(n + 2i - 1) _ /n(n + 1 
(n + 1) 

(n + 1) 
2 

(n + 1) 



+iar 



2A 



i-ll 



tan 



^n + l)A 
2n 



to 



n(n + 1) 



2A 



tan 



i/ Fn+iar 



j-ii- 



n + 1 



2VF, 



n+i ^i-l 



71+1 



(n + 1)A 
2n 



to 



n+i ^i-l 



(100) 



which after substituting for to gives 



l)ar 



i-2 



n(n + 1) 



(n+1) 



(n + 2i-l) Fn+i I 



la 



("-!) 

n+1 



n(n + 1) 



(n+1) 



(n + 1) 2-^F„_|_i y n(n + l)Fn+i — 2A 



(101) 



Equating the power of r on the two sides gives i = 2n/(n — 1) (as in the A = case). Since i should be 
an integer greater than one we find that these conditions are satisfied only for n = 2 (i = 4) and for n = 3 
{^ = 3). 

Consider n = 3; in this case the above equation can be written as 



SJE^a-" + 2 



6 



6F4 - A /F^ 



^ +4F4)a2 + 4 



6 



with the constraint that 



< a < 



6F4 - A 



F&a + 



6 



6F4 - A / F4 



F'i 







I2F4 - 2 A 2F4 



(102) 



(103) 



If we define o = ^/F^b and Fq = F^ \/ 6F4 — A^ the equation can be written in the simplified form 

2&2(46 + ^fli) + (26 + ^6^^ = (104) 

with the requirement that < 6 < — y^3/8^. If we further define Y = —6/^ and r] = —1/^ the equation 
becomes 

2Y^{AY -V6)+ ri{2Y - V^f = 0. (105) 

For a naked singularity to form this equation for Y should have a solution subject to the constraint 
< y < and t] > 0. 

Using the conditions, as mentioned earlier, for the roots of a general cubic we can find the conditions 
on T] for which the above cubic has at least one real root in the desired range. It is found that for 
< ?7 < \/6(-ll + 5^/5)/4 all the thr ee roots are real and at least one of these satisfies < y < \/3/8. 
For 7] > \/6(— 11 + 5\/5)/4 the real root is negative. The range of rj found above implies that for £^ < 
4/^/6(11-5^/5) one gets a naked singularity. This shows that the critical case is also similar to the A = 
case except that the allowed range for ^ has shifted. 
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A similar analysis can be carried out for the case where n = 2 and i = A. By defining a = F^h and 
Fg = F^^^/2>F^ — A one gets a fourth order equation in b. If one subsequently defines Y = —b/^ and 
rj = one gets the equation 

4Y'^{3Y -V3)-r]{Y -Vsf = (106) 

with the consistency conditions < y < l/\/3 and t] > 0. It is found that the above conditions are 
satisfied for < r? < (1590 - 918\/3)/(-9 + 5a/3) « 0.066642 or in terms of conditions on ^ we get 
C < (9 - 5^3)7(1590 - 918^/3) ^ -15.0056. 

3.3 Formation of Trapped Surfaces 

As in the A = case, we now look at the formation of trapped surfaces. Considering the expansion of 
outgoing radial null geodesies as in (ISSp we find that 



nR' / F{r) 2Afl2 

~R~ I ~ Y ~ n{n+l) 



K"-. (107) 



From this it is seen that the condition for trapping, = 0, is met when 

F{r) 2KR^ 



n(n + 1) 

which for n = 2 (4-dimensions) reduces to the well known result 

Fir) Ai?2 



1 (108) 



R 



1. (109) 



It is easy to see that, as in the A = case, for the central shell, trapping coincides with singularity 
formation. 



3.4 Exterior Solution with a negative Cosmological constant 

As before we take the metric in the exterior to be 

ds^ = -f{x)dT^ + g{x)dx^ + x^dQ'^. (110) 

The components of Einstein tensor are the same as in (j39|) -(|41|). Solving the vacuum Einstein equations 
G^iu - ^g^lv = we find 

^ n(n + l)x"-^ 

2A2;"+i + n(n + l)j;"-i + Cn(n + 1) ' ^ ^ 

^ n(n + 1) 

With this the metric in the exterior becomes 

ds'' = - (l + + J^^^\dT^ + (l + + 'dx2 + xW. (113) 

V x'^-i n(n + l)y V X n(n + l)y ^ ^ 

Here C is a constant of integration which is fixed by matching the exterior solution to the interior solution 
at the boundary in exactly the same way as for the A = case and the result is C = —F{rs) where is 
the boundary of the dust cloud. Thus in 4-dimensions the exterior is 

ds^ = -[l + — - ]dT^ +[l + — - dx^ + x^dn^. (114) 
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3.5 The absence of a self-similar solution in the Presence of a A 

It is interesting to note that it is not possible to have a self-similar solution in the presence of a cosmological 
constant. To see this we begin by noting the condition that dimensionless functions made from the metric 
are functions only of t/r continues to hold. This follows because the cosmological constant term in the 
Einstein equations can be absorbed into the energy momentum tensor in the right hand side, by taking 
the A-term as a perfect fluid with equation of state p = —p. Eqn. (j58p then continues to hold, with 
the understanding that the contribution of the cosmological constant is included in the energy-momentum 
tensor. The remaining argument, leading to the conclusion that R is a function of z then follows. 

Now if we have a self-similar solution then we can write R = rR with R being dimensionless. If we 
define j = z then the condition of self-similarity implies that R being dimensionless should be a function 
only of z. With this if we now consider the equation Gn — l^gu = kTh we get 

n(n - 1) 4 



The explicit presence of r in the above equation implies that R cannot be expressed as a function of z 
alone and thus we do not have a self-similar solution in the presence of A. 

The same conclusion also follows from Eqn. ()6ip . With dust matter, the only contribution to pressure 
is coming from the cosmological constant, and this pressure is constant. The Lie derivative on the left hand 
side is thus zero, whereas on the right hand side the presure is non-zero, leading to a contradiction and 
showing that such a Kiling vector field cannot exist. Physically speaking, the presence of a cosmological 
constant introduces a length scale which prevents self-similarity. 

4 Conclusions 

We have studied the collapse of inhomogeneous spherically symmetric dust distribution in arbitrary number 
of space dimensions both in the absence and in the presence of a cosmological constant. From the analysis 
presented we see that even though naked singularity is allowed in all dimensions there is more freedom on 
initial conditions for obtaining naked singularity in 2-1-1, 3-1-1 and 4-1-1 dimensions, both in the absence as 
well as the presence of a negative cosmological constant. We have also seen that the formation of trapped 
surfaces is similar in all dimensions with the central shell getting trapped at the same time when it becomes 
singular. For outer shells trapping occurs before those shells become singular. We also saw explicitly that in 
the absence of a cosmological constant, globally naked self-similar models can be constructed all dimensions, 
whereas in the presence of a cosmological constant such a solution cannot be constructed. 

In the second paper in this series, we will study quantum field theory on the curved background provided 
by the classical solutions presented here, including the emission of Hawking radiation from an n-dim AdS 
black hole. In a third paper we will carry out a canonical quantization of this model, and also address the 
issue of black hole entropy, following the methods of [37j . 
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